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ABSTRACT 

New Vic CCD photometry, obtained with integration times of 20s, of the 
sdOB+degenerate-dwarf eclipsing binary system AA Dor has provided new complete 
light curves with an rms scatter about a mean curve of ±0.004 mag. These data are 
analysed with an improved Light2 light curve synthesis code to yield more accurate 
determinations of the radii of both stars, the orbital inclination, and the flux ratio 
between the two components. These radii are only a little different from the values 
' derived 25 years ago from less complete data, but the uncertainties on these values are 

, improved by a factor of two. The apparent discrepancy remains between the surface 

gravity of the sdOB primary star obtained from the light-curve solution with the pub- 
lished spectroscopic orbit and that obtained from NLTE analysis of high-resolution 
spectra of the sdOB star. 

The substantial reflection effect in the system is adequately represented by the 
Jjj | Light2 code with a bolometric albedo of unity in light curves extending from 0.35 \x m 

to 2.2 /im. However, there arc differences at individual wavelengths in the derived 
albedo, which may indicate redistribution of flux from shorter wavelengths into the V 
and Ic passbands. 

•rH . 

. Key words: stars: binaryeclipsing; stars: binaryclose; reflection effect; stars: astro- 

' physical parameters; stars: individual (AA Dor). 
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1 INTRODUCTION 

The blue star LB 3459 (HD269696; V ~ 11.1m, RA = 
05' l 31 m 40 s .3, Dec = -69° 53' 02" .2 (2000)), and situated 
at high galactic latitude (b = —32°. 18), was discovered 
to be variable in brightness by Kilkenny, Hilditch & Pen- 
fold (1978). Observations with the University of Cape Town 
high-speed photometer revealed the light curve of an eclips- 
ing binary (designated AA Doradus) with an orbital period 
of 0.261 day and a primary eclipse of depth 0.4 mag in the 
V-band. Spectroscopic and further photometric observations 
followed over many years, and the present understanding of 
the properties of this remarkable binary system is summa- 
rized in Hilditch, Harries & Hill (1996) (hereafter HHH), and 
in Rauch (2000) where references to previous publications on 
this system may be found. 

HHH obtained echelle spectra of AA Dor from the 
Anglo- Australian Telescope and determined a velocity semi- 
amplitude for the sdOB primary star of 40.8 ± 0.7 kms~ . 
Many attempts to detect spectral signatures of the sec- 



ondary component from these data failed, and HHH were 
forced to adopt a representative mass of 0.5 Mq for the sdO 
star in order to calculate the mass of the secondary star from 
the derived mass function and the earlier light-curve analy- 
ses, and establish the absolute sizes of the stars in this sys- 
tem. The sdOB star seemed to be an unremarkable member 
of that class of evolved stars, whilst the unseen secondary 
component at 0.09 Mq seemed to be a degenerate red dwarf 
of surface temperature ~ 2000A'. The derived surface grav- 
ity of the primary star was found to be log g = 5.53 ± 0.03 
which was in some disagreement with an earlier NLTE anal- 
ysis by Kudritski,e£ al. (1982), where a value of 5.3 ± 0.2 
was derived. HHH recommended that a new determination 
of log g from spectral line profiles would help to resolve this 
issue and permit a direct determination of the mass of the 
sdOB star and therefore also the secondary component. 

The more recent work by Rauch (2000) addressed this 
issue directly by detailed NLTE analysis of high resolution 
optical and ultraviolet spectra. Rauch obtained an effective 
temperature for the primary component of 42000 ± 1000 K, 
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in agreement with Kudritski et al., and a value of log g = 
5.21 ± 0.10 which clearly exacerbated the problem. He also 
derived a mass for the sdO star of 0.330± 0.006 from compar- 
ison of the location of the primary in the T e f f vs log g plane 
with evolution tracks for post-RGB stars. The consequent 
mass for the secondary component was 0.066 ± 0.001 Mq 
which places it firmly in the brown dwarf category. 

With the discrepancy about the surface gravity still un- 
resolved, we decided to secure new complete light curves of 
AA Dor in the V and I bands in order to re-determine the 
radii of both stars relative to the semi-major axis of the 
relative orbit. Most recently, Rauch & Werner (2003) have 
published a new spectroscopic study of AA Dor based upon 
many 180s-exposure UVES spectrograms obtained with the 
ESO Very Large Telescope. 



2 PHOTOMETRIC OBSERVATIONS 

Because the primary eclipse of AA Dor has a total dura- 
tion of only 26 minutes and a depth of 0.4 mag, it is clear 
that photometry has to be secured with excellent time res- 
olution in order to ensure that the detailed shape of the 
eclipse curve is established without compromise. Accord- 
ingly, the University of Capetown CCD camera, mounted 
on the 1.0-m telescope at the South African Astronomical 
Observatory, was used to secure these data. This camera op- 
erates in frame-transfer mode, so that there is essentially no 
time lost between exposures. Integration times of 15- or 20-s 
were used for the V observations, and 20-s for the Ic obser- 
vations, ensuring that the eclipse curves were well sampled. 
To be certain of this time resolution, all the observations on 
a single night were obtained in one filter only, rather than 
alternating between the V and Ic passbands. The field of 
view of the UCT CCD camera on the 1-m telescope is only 
about 1.8 x 1.2 arcminutes so that, despite A A Dor being a 
foreground object to the Large Magellanic Cloud, the field is 
remarkably free of useful comparison stars. One comparison 
star is ~ lm fainter than AA Dor at maximum, and the 
rest are all at least 2m fainter. 

Observations of the AA Dor field were obtained on 4 
consecutive nights, two nights with incomplete light curves, 
and two nights with complete light curves, one in each of the 
V and Ic passbands that are in routine use at the SAAO. 
The data were reduced with an SAAO version of DoPhot, 
slightly modified to reduce data "on-line" so that an al- 
most immediate assessment of photometric quality can be 
made. Differential magnitudes were derived relative to the 
brightest comparison star using standard SAAO software. 
The constancy of this star was established by comparing it 
with the fainter stars in the field. 

Infrared JHK photometry was obtained using the 
Mklll photometer (Glass 1973), with the //50 chopping 
secondary, on the 1.88-m reflector of the SAAO. Observa- 
tions were made during two allocations of time: 1984 Jan- 
uary 17-23 and 1984 January 27-30. For each filter, obser- 
vations were obtained in modules of two 10-s integrations. 
Beamswitching occurred twice in each module, which was re- 
peated as necessary in order to achieve the desired accuracy. 
On a given night, a single filter was selected and A A Dor 
monitored for an entire orbital period. The star HD37027 
was selected as a comparison star since it was an early-type 



star close to A A Dor and which Cousins (1987) subsequently 
published as a secondary UBV standard star. 

Infrared standard stars were selected from the SAAO 
list, which Carter (1993,1995) later published, and observed 
in so far as was practicable at approximately the same air- 
mass as AA Dor. These enabled differentially- corrected in- 
frared magnitudes to be placed on the standard system. The 
intention was to detect a signature of the cool secondary 
component in the infrared light curves, specifically by ex- 
amining whether the depth of the secondary eclipse would 
be found to dip below the occultation of the reflection effect 
by the primary, that is, below the level recorded immediately 
outside the primary eclipse, and thereby yield some measure- 
ment of the infrared flux from the secondary star. The data 
were of poorer quality than had been hoped, and they failed 
to show a clear signature of the flux from the secondary 
star beyond the known reflection effect; they are included 
graphically as a further test of reflection-effect modelling, 
considered in section 3.3. 



3 ANALYSIS OF THE LIGHT CURVES 

3.1 The Vic data 

The most recent orbital ephemeris for AA Dor has been 
determined by Kilkenny et al. (2000) from times of pri- 
mary eclipse minima observed from 1977 to 1999. This lin- 
ear ephemeris has a reference time of primary minimum 
of HJD 2443196.34868 ± 0.00002, and an orbital period of 
0.261539731 ± 0.000000002 d, and has been used to con- 
vert the heliocentric Julian Dates into orbital phases. To 
improve photometric accuracy, the individual observations 
were added in orbital-phase bins of A cf> — 0.0020 and then 
differential magnitudes formed between the variable and the 
comparison star. Most bins had 2-3 observations yielding a 
typical Poisson uncertainty of ±0.004 m in V and ±0.003 m 
in Ic- 

The two complete light curves are plotted in Figure 
with a total of 500 binned observations defining the V curve 
from 1157 individual observations obtained on JD2451880, 
and in Figure[5]with a total of 500 binned observations defin- 
ing the Ic curve from 1287 individual observations obtained 
on JD2451879. On each of these nights complete light curves 
were secured and the small region of orbital-phase overlap 
with data obtained at the beginning and end of each night 
showed no differences. The overall V curve shows some mi- 
nor irregularities from a symmetric curve, but the Ic curve 
is entirely symmetric. 

3.2 Geometric parameters 

The two stellar radii, expressed in terms of the semi-major 
axis of the relative orbit, and the orbital inclination are de- 
termined very well from the annular primary eclipse curve 
alone, as noted in HHH. We have used the LIGHT2 light- 
curve-synthesis code (Hill 1979, Hill & Rucinski 1993) to 
analyse the complete V and Ic curves separately, and the 
results are given in Table We adopted standard values 
of linear limb darkening for both stars, together with stan- 
dard values of the gravity darkening exponent. A mass ra- 
tio of 0.20 (secondary/primary) was adopted from Rauch 
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Figure 1. The complete V band light curve of AA Dor defined 
by 500 binned observations (crosses with Poisson uncertainties of 
±0.004 m). The solid line shows the LIGHT2 model fit to these 
data, whilst the dotted line is the theoretical curve with the 
albedo fixed at unity. The observed-calculated (O — C) values 
are also plotted displaced by +0.45 m from zero. 
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Figure 2. The complete Iq band light curve of AA Dor defined 
by 500 binned observations (crosses with Poisson uncertainties 
of ±0.003 mag). The solid line shows the LIGHT2 model fit to 
these data, whilst the dotted line is the theoretical curve with 
the albedo fixed at unity. The observed-calculated (O — C) values 
are also plotted displaced by ±0.50 m from zero. 

(2000). The effective temperature of the sdO star was taken 
to be 42000J\" also from Rauch's (2000) NLTE spectral anal- 
ysis, and that of the cool secondary star to be a nomi- 
nal 2000-K\ Solutions of the Zc-band curve with different 
fixed values of the secondary temperature (T sec ) showed 
that the \ 2 value increases by a factor 2 over the range 
2000if < T sec < 6000T\\ The reason for this lack of depen- 
dence on the secondary temperature is simply that the flux 
ratio (primary /secondary) is 10 4 in the Ic band and 10 5 
in the V band immediately after egress from the primary 
eclipse, when we see only the averted unheated hemisphere 
of the secondary. At the maximum of the reflection effect 



Table 1. Geometric parameters 



passband 




^sec 


incl 


rms 


V 


0.1413 


0.0759 


89.60 


±0.0053 




±0.0008 


±0.0006 


±0.56 




Ic 


0.1425 


0.0769 


88.82 


±0.0039 




±0.0007 


±0.0005 


±0.20 




average 


0.1419 


0.0764 


89.21 
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Figure 3. The V band primary eclipse light curve of AA Dor 
(crosses with Poisson uncertainties of ±0.004 m). The solid line 
is the LIGHT2 model fit to these data solving for the geometrical 
parameters and the albedo of the secondary star. The dotted line 
is a theoretical curve calculated from the average geometrical pa- 
rameters and adopting an albedo of unity for the secondary.Thc 
observed-calculated (O — C) values are also plotted displaced by 
±0.45 m from zero. 



this ratio has reduced to 13 in Ic and 19 in V, which sug- 
gests that it should be possible to detect the heated face 
of the secondary in near-infrared spectra obtained at the 
appropriate orbital phases. 

It is clear from the results in Table 0that the solutions 
for the mean, or volume, radii and the inclination from the 
V and Ic curves are in excellent agreement, and that the 
Ic solution is slightly more precise. The average values have 
been adopted in all subsequent analyses as fixed parameters, 
and we should note that they are little different from those 
found in earlier studies from less complete and less accurate 
data, but they are more precisely determined. 

The excellent agreement between these V and 7c-band 
observations and the model light curves are shown in detail 
through the primary and secondary eclipse curves in Fig- 
ure [|J Figure |1J Figure and Figure |S| V and Ic curves 
defined by the individual observations were also solved, and 
the same values for the geometrical parameters were deter- 
mined as in Table 1 but with larger uncertainties. Combin- 
ing the observations into these phase-binned values has not 
systematically affected the derivation of the relative stellar 
radii. 
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Figure 4. The Iq band primary eclipse light curve of AA Dor 
(crosses with Poisson uncertainties of ±0.003 m). The solid line 
is the LIGHT2 model fit to these data solving for the geometrical 
parameters and the albedo of the secondary star. The dotted line 
is a theoretical curve calculated from the average geometrical pa- 
rameters and adopting an albedo of unity for the secondary.The 
observed-calculated (O — C) values are also plotted displaced by 
+0.50 to from zero. 
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Figure 6. The Iq band secondary eclipse light curve of AA Dor 
(crosses with Poisson uncertainties of ±0.003 to). The solid line 
is the LIGHT2 model fit to these data solving for the geometrical 
parameters and the albedo of the secondary star. The dotted line 
is a theoretical curve calculated from the average geometrical pa- 
rameters and adopting an albedo of unity for the secondary.The 
observed-calculated (O — C) values are also plotted displaced by 
+0.050 m from zero. 
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Figure 5. The V band secondary eclipse light curve of AA Dor 
(crosses with Poisson uncertainties of ±0.004 m). The solid line 
is the LIGHT2 model fit to these data solving for the geometrical 
parameters and the albedo of the secondary star. The dotted line 
is a theoretical curve calculated from the average geometrical pa- 
rameters and adopting an albedo of unity for the secondary.The 
observed-calculated (O — C) values are also plotted displaced by 
+0.050 m from zero. 



3.3 Reflection effect 

The manner in which the reflection effect is modelled does 
not significantly affect the solution for the geometrical pa- 
rameters in A A Dor, which is why that part of the analy- 
sis has been separated. In HHH, it was explained that the 
LIGHT2 code had difficulty in fitting the observed reflection 
effect, and recourse was made to rather extreme measures 
to fit the data, by changing the limb-darkening coefficient 
of the heated secondary star to negative (limb-brightening) 
values. Whilst such drastic action did provide a solution, it 
was clearly not satisfactory. 



With these new complete light curves, it was necessary 
to resolve this issue, particularly in view of the fact that 
users of the Wilson-Devinney (wd) code (Wilson & Devin- 
ney 1971, Wilson 1979, 1990) on eclipsing binaries with simi- 
lar large-amplitude reflection effects did not experience such 
major difficulties (e.g. Wood et al. 1993, Zola 2000), al- 
though it did seem necessary for Wood et al. to adopt some 
unusually low values for the limb-darkening coefficients. A 
detailed reexamination of the code by GH finally revealed a 
one-line error in the reflection calculation subroutine, a miss- 
ing factor 7r coupled with an incorrect expression involving 
the angle to the surface normal of emergent intensity. Mak- 
ing this correction increased the calculated amplitude of the 
reflection effect up to that calculated by the WD code. 

This revised code was tested against a set of light curves 
kindly provided by Professor R.E.Wilson (private communi- 
cation). These monochromatic light curves were generated 
with the current version of the WD code using nominal values 
of the two radii, the orbital inclination, and the tempera- 
tures of the two stars in AA Dor, and setting the bolometric 
albedo of the secondary star to be unity. The monochro- 
matic wavelengths were set to the mean wavelengths of 
the Stromgren-Crawford uvby filters, the V and Ic filters, 
and the infrared JHK passbands. Since it is already long- 
established that the WD and LIGHT2 codes agree in all geo- 
metrical respects, (cf., for example, Bell et al. 1991), these 
uvby and Vic light curves were solved individually by the 
revised LIGHT2 code for the albedo of the secondary compo- 
nent. Two descriptions of the reflection effect were used, one 
adopting that in the WD code of complete absorption and re- 
radiation of the incoming flux from the primary star by the 
secondary star (i.e. heating), and the other that provided by 
Sobieski(1965) where some of the incoming flux may be scat- 
tered by free electrons in the atmosphere of the secondary 
and the remainder goes into direct heating. The direct heat- 
ing approach yielded a mean albedo from the 9 light curves 
of 0.97 ± 0.06, whilst the Sobieski approach yielded a mean 



Table 2. Values of heating albedo for the secondary component 



passband heating — only Sobieski 



u 


0.46 


0.49 




±0.12 


±0.16 


V 


0.45 


0.48 




±0.06 


±0.06 


b 


0.39 


0.42 




±0.05 


±0.06 


V 


0.55 


0.63 




±0.07 


±0.10 


V 


1.43 


1.70 




±0.05 


±0.06 


Ic 


2.02 


2.40 




±0.04 


±0.05 



albedo of 1.16 ± 0.07. This direct comparison between the 
WD and the LIGHT2 code shows a satisfactory agreement in 
the representation of the reflection effect. 

The V and Ic light curves were analysed with the 
LIGHT2 code, adopting black-body fluxes convolved with the 
standard filter response curves, and solving for the two radii, 
the orbital inclination, and the albedo (a sec ) of the sec- 
ondary component. Because of the large range in temper- 
ature across the heated hemisphere of the secondary com- 
ponent (2,000 — 20, 000 AT), we used the option in light2 
to use a limb-darkening coefficient appropriate to the lo- 
cal temperature of each code pixel on the secondary surface 
rather than a mean value. Both modes of calculating the re- 
flection effect were used, and these provided identical values 
of the geometrical parameters. These values are reported in 
Table as already noted. The values for a sec are given in 
Table|5]for the two wavelengths and both modes of calculat- 
ing the reflection effect. The root-mean-square (rms) scat- 
ter of the binned observations about the fitted light curves 
(noted in Table are only 0.001 m larger than the Poisson 
uncertainties on these binned data, demonstrating that the 
fits are very good. 

The uvby light curves of AA Dor published by Kilkenny, 
Penfold & Hilditch (1979) (Figure[7J were also analysed with 
the LIGHT2 code solving only for the albedo of the secondary 
component in both the heating-only mode and in the So- 
bieski mode, and adopting as fixed values the average ge- 
ometrical parameters given in Table 1 from the V and Ic 
light curve solutions. These results, together with those de- 
termined from the V and Ic curve solutions are listed in 
Table H 

Irrespective of the mode of reflection calculation used, 
there are significant departures from an idealised a sec = 1.0 
model. In Figures 3-7, these solution curves are plotted 
as solid lines, there being no difference between the two 
modes. Also plotted in these figures are theoretical light 
curves for a binary with the average geometrical parame- 
ters and a 3e c = 1.0. The uvby light curves give an aver- 
age value of a sec — 0.50 ± 0.05 for the Sobieski mode and 
otsec = 0.46±0.04 for the heating-only mode. The a sec — 1.0 
light curves systematically do have reflection effects with 
larger amplitudes than the observations, but those differ- 
ences amount only to ~ 0.01 — 0.02 m, and the observational 
uncertainties are of the same order per observation. For the 
better-defined V and Ic curves, the observations present a 
reflection effect of higher amplitude than the a sec = 1.0 
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Figure 7. The uvby light curves of AA Dor (crosses) published 
by Kilkenny et a/.(1979), and the LIGHT2 model fit to these data 
(solid lines) adopting the average geometrical parameters and 
solving for the secondary albedo. The dotted lines show the the- 
oretical curves for an albedo of unity. 

curves and these systematic differences are clearly evident, 
particularly in Figure [S] and Figure |S| Since these single 
epoch V and Ic light curves were obtained 22 years after 
the uvby data, we cannot rule out the possibility of some 
intrinsic variability in the system that conspires to simulate 
an enhanced reflection effect in two wavelengths and a sup- 
pressed effect at shorter wavelengths. Alternatively, perhaps 
this is observational evidence for some redistribution of flux 
from the near-UV and blue regions into the visual and near- 
IR regions by means of substantial line blocking/blanketing. 
Note that this intrinsically very cool secondary star has an 
unperturbed T e ff = 2000 K, and a temperature at the sub- 
stellar point directly underneath the sdOB primary star of 
some 20, 000 if. It would be interesting to know what type 
of overall spectrum would be calculated for such an illumi- 
nated hemisphere, (c/., e.g., Barman, Hauschildt & Allard, 
2002). In our analysis, the black-body fluxes are a reasonable 
representation of the broad-band sdOB spectrum. But the 
spectrum of the heated hemisphere of the secondary com- 
ponent is likely dominated by emission lines and our black- 
body representation is clearly a limitation. It might also be 
expected that circulation currents would exist on the facing 
hemisphere of the irradiated secondary star (Podsiadlowski 
1991), and that the limb-darkening values would be reduced 
compared to normal. Reducing the limb-darkening values in- 
creases the calculated amplitude of the reflection effect for a 
given value of the albedo, as already demonstrated by Wood 
et al, and by HHH. 

In Figure |HJ the infrared JHK band data are plotted, 
together with theoretical curves calculated from the average 
geometrical parameters and with a sec = 1-0. Despite the ob- 
vious large uncertainties in these data, such theoretical light 
curves are consistent with these data and do not suggest any 
large-scale departure from a model with a bolometric albedo 
of unity for the secondary star. 

We note again that Wood et al. experienced some diffi- 
culties in fitting their multi-colour observations of HW Vir 
with the WD code in the sense that they had to adopt low 
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Figure 8. The JHK light curves of AA Dor (crosses) obtained 
by AELG. The solid lines show the theoretical curves calculated 
from the average geometrical parameters for an albedo of unity. 

values of the limb-darkening coefficients in order to increase 
the theoretical amplitude of the reflection effect up to that 
observed. In a study of the system PG1336-018, Zola (2000) 
modified the WD code to assign limb-darkening coefficients 
appropriate to the local temperatures of each of the stel- 
lar surface pixels in the code, and obtained an albedo value 
in the range 0.45 — 0.95 dependent on other parameters for 
the system. Thus there is evidence that in these binary sys- 
tems displaying extreme reflection effects the accepted stan- 
dard light-curve-synthesis codes experience some difficulties 
in matching the observations. 



4 ASTROPHYSICAL PARAMETERS 

The results of several observational programmes investigat- 
ing the fraction of hot subdwarf stars that are found in bi- 
nary systems have been published in recent years. Specifi- 
cally, Maxted, Heber, Marsh & North (2001) and many ref- 
erences therein, have demonstrated that the observed frac- 
tion of extreme horizontal branch stars that are members 
of binary systems with orbital periods less than 10 days is 
~ 70%. Two recent papers by Han et al. (2002) and Han, 
Podsiadlowski, Maxted & Marsh (2003) have discussed the 
various evolution channels that can explain the observed 
properties of sdB,sdOB, and sdO stars in binary systems, 
concluding that the different channels can explain why some 
hot subdwarfs are found in short-period binaries, others in 
binaries with much larger orbits, and some are now merged 
single objects. The first common-envelope-ejection (CEE) 
channel is the evolution path where the initially more mas- 
sive star experiences dynamical mass transfer on its first 
ascent of the red-giant branch, leading to the formation of 
a common-envelope followed by a spiral-in phase with the 
binary system attaining a short orbital period (P < Wd) 
and being composed of a hot subdwarf and a main-sequence 
companion. The distribution of masses of the hot subdwarf 
stars in this evolution channel is expected to lie in the range 
0.33 — 0.47 Mq, with a main peak at around 0.46 Mq and a 
much smaller peak around 0.33 Mq. These binary popula- 



Table 3. Astrophysical parameters 



parameter primary secondary 



mass Mq 


0.33 


- 0.47 


0.064 - 


0.082 


radius Rq 


0.179 


- 0.200 


0.097 - 


0.108 


log g (cgs) 


5.45 


- 5.51 


5.27 - 


5.29 


V r ot km s _1 


34.7 


- 38.6 


18.7 - 


20.8 



tion synthesis models do succeed in explaining the observed 
distributions of properties of sdB,sdOB,sdO star binaries. 
Specifically with regard to AA Dor, the first CEE models 
seem to provide a sensible explanation for its observed prop- 
erties. 

Since the spectrum of the secondary star in AA Dor 
has not yet been observed, so that we do not know its or- 
bital velocity curve, we have to assume a mass for the sdOB 
star, and then calculate all the absolute astrophysical pa- 
rameters from the derived values of relative radii, orbital 
inclination, and the velocity semi-amplitude of the primary 
component. HHH derived a value of K pr i — 40.8±0.7 km s _1 
whilst Rauch & Werner (2003) determined a value of K pr i = 
39.19 ± 0.05 fcms" 1 . Adopting an expected range of subd- 
warf mass of 0.33 — 0.47 Mq, with both values of K pr i, and 
the geometric parameters from Tabled we obtain the ranges 
of possible values of masses, radii, surface gravities, and syn- 
chronous rotational velocities given in Table [3] 

The corresponding range of mass ratio is 0.168 — 0.200 
and has no effect upon the derived relative radii, orbital 
inclination etc from the light curve analyses. These ranges 
result from the assumed range of primary star mass, and 
are larger than the range of uncertainties that would result 
from propagation of the formal uncertainties of the derived 
quantities. They may be compared with the values deter- 
mined by Rauch (2000) and Rauch & Werner (2003) from 
their analyses of spectra of AA Dor. 

Rauch (2000) derived a mass for the sdOB star of 
0.330 ± 0.006 Mq from comparison of their derived effective 
temperature (T e f / = 42000 ± 1000 K) , and surface gravity 
(log g = 5.21 ± 0.10) with theoretical models for the evo- 
lution of subdwarf stars of different masses plotted in the 
T e f f vs log g plane. This value of surface gravity disagrees 
with that determined from the velocity curve of the primary 
component and the above light curve analysis. We note that 
Rauch & Werner (2003) state: "Since the analysis of Rauch 
(2000) was hampered by the relatively long exposure times 
(lh and 2-3h, respectively) and hence, a relatively large or- 
bital velocity coverage (the observed line profiles were broad- 
ened by the star's rotation as well as by smearing due to 
orbital motion within the observations), it has been specu- 
lated that log g from Rauch (2000) is somewhat too low." 
The new data reported in Rauch & Werner (2003) (here- 
after RW) have exposure times of just 180s and remove this 
problem. From NLTE model atmospheres, RW fitted theo- 
retical line profiles for the H el IX 4686 absorption line and 
derived T e ff = 44000 AT, log g — 5.4, without published un- 
certainties, and Vrot = 43 ± 5fcms _1 . They note, however, 
that this determination of T E ff is not as definitive as that 
determined by Rauch (2000) from the ionization equilibria 
of many species, whilst the value of log g = 5.21 was deter- 
mined from the broad Balmer lines of hydrogen. Later, RW 
note that: " Unfortunately, significant effects of changes in 
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g are detectable only in the outer line wings where the data 
reduction of the broad Balmer lines in the echelle spectra is 
not very accurate." The value of log g = 5.4 for ffe//A4686 
is remarkably close to the range in Table [3] derived from 
our radial-velocity and light-curve analyses, and the value 
of Vrot — 43 ± 5 km s _1 is consistent with synchronism. 

These new results from RW suggest that the apparent 
discrepancy in surface gravity may already be resolved (cf. 
the above results from the H el IX 4686 line), although RW 
did reject this solution in favour of the Rauch (2000) value. 
With log g = 5.21 and T eff = 42000/sT fixed, RW determined 
a rotational velocity of V ro t = 47 ± 5 km s _1 from the A 4686 
line, which is 20 — 35% faster than synchronism. 

The ranges of mass, radius and surface gravity for the 
sdOB star lie within the expectations of the binary popu- 
lation synthesis models of Han et al., whilst the ranges of 
mass, radius and surface gravity of the unseen secondary 
component are consistent with models for the lowest mass 
dwarfs (Dorman et al. 1989) at the boundary of the lower 
main sequence and the substellar brown-dwarf region. 

The referee has commented that the T e / / , log g values 
for the primary component of AA Dor place it in an unusual 
and relatively-short-lived phase of evolution compared to 
most hot subdwarf stars. When compared to the binary- 
star evolution models of Han et al. (2002), or the single-star 
evolution models of Dorman, Rood & O'Connell (1993), its 
location would suggest that the AA Dor primary is evolved 
beyond the terminal-age extreme horizontal branch (taehb) 
towards the white dwarf region, a phase of evolution that 
lasts ~ 10 7 yr. We may ask whether the relative numbers 
of sdO and sdOB stars to all hot subdwarf stars would be 
consistent with that view, at least in the sense of being a 
small proportion of the total. 

From the analysis of Zone 1 of the Edinburgh-Cape 
(EC) Blue Object Survey, Kilkenny et al (1997) found 27 
sdOB stars (like the A A Dor primary), and 20 sdO stars 
out of a total of 357 identified hot subdwarfs, making pro- 
portions of ~ 8% and ~ 6% respectively. The vast majority 
(~ 76%) are sdB stars, with the next most common being 
the He-sdO and He-sdB stars at ~ 11%. From the Ham- 
burg Schmidt Survey (HSS), Lemke et al (1997) selected 400 
stars for more detailed spectroscopic analysis. They identi- 
fied 34 He-sdO stars and 13 normal sdO stars, the latter be- 
ing ~ 3% of the total. Their derived values of T e f j for these 
normal sdO stars lay in the wide range 32, 0007f — 100, 000K , 
whilst the values of log g were in the range 5.0 — 7.0 with 
the average at 5.6 ± 0.6, very similar to the AA Dor pri- 
mary. These statistical results, and the model-atmosphere 
analyses, are consistent with the evolutionary expectations 
that sdO, sdOB stars do lie beyond the taehb in this short- 
lived phase of evolution. By contrast, the very recent spec- 
tral analyses by Edelmann et al. (2003) of sdB stars discov- 
ered in the HSS is in agreement with the work of Maxted et 
al. (2001) on the locations of the sdB stars mostly within the 
bounds between the zero-age and the terminal-age extreme 
horizontal branch. 

5 SUMMARY 

Revised values of the mean radii of the two stars in AA Dor 
have been determined from these new Vic light curves 



which differ only a little from those determined earlier, but 
they are defined more precisely. These results together with 
the published orbital parameters from HHH and RW provide 
new determinations of the astrophysical parameters with an 
adopted range of primary star mass. These results are com- 
patible with the theoretical evolution models for binary sys- 
tems containing hot subdwarf stars published by Han et al. 
(2002, 2003). 

An error in the reflection effect subroutine within the 
LIGHT2 code has been corrected with the result that it agrees 
with the accepted standard WD code. Analysis of the ob- 
served reflection effect demonstrates significant departures 
within some passbands from an idealised albedo of unity for 
the heated hemisphere of the secondary star, but the overall 
bolometric value is most likely close to unity. It is suggested 
that further theoretical work should be done on modelling 
the reflection effect in these types of binaries containing stars 
with very different effective temperatures. 
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